Hawking Radiation from an Acoustic Black Hole on an Ion Ring 
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In this Letter we propose to simulate acoustic black holes with ions in rings. If the ions are 
rotating with a stationary and inhomogeneous velocity profile, regions can appear where the ion 
velocity exceeds the group velocity of the phonons. In these regions phonons are trapped like light 
in black holes, even though we have a discrete field theory and a nonlinear dispersion relation. We 
study the appearance of Hawking radiation in this setup and propose a scheme to detect it. 
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In 1974 Hawking showed that the theory of quan- 
tum fields in curved spacetime predicts that, surpris- 
ingly, black holes emit thermal radiation [1]. Unfortu- 
nately, the temperature of this radiation is too small to 
be detected for typical astrophysical black holes. Fur- 
thermore, the original theoretical derivation suffers from 
the problem that the wave equation is assumed to be 
valid on all scales, whereas the theory of quantum fields 
in curved space is assumed to break down at the Planck 
energy. Unruh showed that the Hawking effect is also 
manifested in analogous hydrodynamical systems which 
have a region of supersonic flow and hence a sonic horizon 
0. Such analogous systems offer great advantages, since 
the effect can potentially be accessible to experiments. 
Moreover, its robustness can be examined based on the 
well known microphysics of the hydrodynamical systems. 
This will contribute to deepen our understanding of the 
Hawking effect also in gravitational black holes. 

The hydrodynamic analogy of gravitational spacetimes 
has inspired many proposals for experimental tests of 
Hawking radiation in continuous fields in recentyears Q , 
e.g., phonons in Bose-Einstein condensates y-roj , Fermi 
gases 0, superfiuid Helium J^, slow light |9lllo|, and 
nonlinear electromagnetic waveguides [ll|. So far, no 
proposal has been physically implemented. 

In the present work we show how to build an analog 
model of a black hole in an experimentally realizable sys- 
tem of ions. This is the first experimental proposal in 
which a discrete system is completely analyzed in the 
discrete limit (see also A sublinear dispersion re- 

lation at high wave numbers naturally results from the 
discreteness of the physical system. This affects the tra- 
jectories of blue-shifted waves close to the event horizon 
1^ . The dispersion relation is, additionally, nontrivial 
at low wave numbers because of the long range Coulomb 
force. We study how much this affects the appearing 
Hawking radiation. Explicit numerical calculations show 
that the Hawking effect is robust against such short scale 
modifications, e.g., for a continuous field with a sublinear 



dispersion relation [iSj and a discretized field on a falling 
lattice [13] . Our proposal uses a parameter regime which 
is accessible in experiments at temperatures currently 
achieved. Thus, it could lead to the first experimental 
observation of Hawking radiation. 

The main idea of our proposal can be summarized as 
follows. We are constructing a discrete analog of a hydro- 
dynamical system with super- and subsonic regions in a 
quadrupole ring trap 1^ as schematically depicted in 
Fig. [1] Ions rotate on a ring with circumference L with 
an inhomogeneous, but stationary velocity profile v{9). 
Thus the ions are inhomogeneously spaced. The neces- 
sary (de-)acceleration of the ions is guaranteed by addi- 
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FIG. 1: Schematic depiction of the ion rotation with velocity 
V (9) and the phononic excitations with velocity ±c{6), that 
depend on the varying ion spacing. The external force {0) 
localized at the transition between the super- and subsonic 
region (de-) accelerates the ions. The inset shows a typical 
velocity profile with a black hole like region for v (6) > c (9). 



2 



tional electrodes exerting a force F'^ (9) on the ions. The 
harmonic oscillations of the ions around this equilibrium 
motion are phonons with velocities c{9) (x {v (9)) 
Regions with large enough ion spacings and low enough 
velocities are supersonic; here phonons can only move 
into the direction of the ion flow and are trapped like 
light inside a black hole. We consider a system with a 
super- and a subsonic region. The border between these 
regions is analogous to a black hole horizon and will be 
shown to emit Hawking radiation. 

In the following we explain the detailed setup of our 
proposal. The dynamics of N ions with mass m and 
charge e are described by the Hamiltonian 
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with the Coulomb potential and a local external po- 
tential [9] . Instead of specifying we will impose 
an angular velocity profile v (9) by fixing the equilibrium 
positions. The required is then determined through 
the difference between the ion acceleration 9^ {t) and the 
Coulomb force i^f (6/?, . . . , 9%) . 

We impose the classical equilibrium positions 



9Ut) = g{t/N + t/T), 



(2) 



where g maps the normalized indices i/N G [0, 1] mono- 
tonically increasing onto the angles 9 ^ [0, 2tt] and is 
periodically continued. The angular velocity profile is 
v{9) = g' (5-1 (9)) /T. We make the choice 
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and h[s) 



with a — and p = 2 

15/8s — 5/4:3^ + 3/8s^. v (9) — v„iax is constant in the 
supersonic region aVminT ^ 6* ^ 1 — crVminT and v (9) = 
Vmin constant outside of it {vmin < Stt/T < Vmax)- We 
choose Vmin = 27r • 0.83/T and that the small transition 
regions 271,272 contain 0.04A^ and O.IN ions. 

We treat small perturbations around the equilibrium 
motion 9i (t) = 9^ (t) + 59i (t) and expand the Hamilto- 
nian to second order in 69 j 
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with the time dependent force matrix fij and the canon- 
ical operators S9i and —ihdsg^ describing the phononic 
oscillations of the ions. The quasifree quantum dynam- 



ics of this harmonic system are governed by the classi- 
cal linear equations of motions for the first and second 
moments. We have checked that anharmonic effects do 
not significantly alter the results presented in this Letter. 
With Floquet analysis we find systems with stable and 
unstable dynamics in the interesting parameter regime 
similar to In the presence of two horizons, black 
hole laser instabilities ITj may contribute to this. How- 
ever, the instabilities become dominant at such late times 
that they are not important for the simulations presented 
here. Nevertheless, the examples given in this Letter rep- 
resent stable systems. 

In order to get some insight, we consider the limit of 
an infinite number of ions, where we can make an anal- 
ogy with the standard Hawking effect for the scalar field 
$ (6»f' {t) , t) = 59.1 (t) as observed in Then for a slowly 
varying v (9) and Coulomb interactions between neigh- 
boring ions only the system Lagrangian becomes 
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with the dispersion relation D {9, k) = c{9) k + O (fc'^) , 
the conformal factor p{9) ^ n{9) ■ mL'^ / {2tt)'^ , and the 
density n{9) = N/{v{9)T). This scalar field satisfying 
a linear dispersion relation at low wave numbers with 
sound velocity c (9) = y^2{2TT)^n (9) e^/ (47reo • mL-^) is 
analogous to a massless scalar field in a black hole space- 
time [2]. Its quanta cannot escape a supersonic region 
with V (9) > c (9) like photons trapped inside a black 
hole. The horizon of this analog model is located at 
c(9h) = v{9h) with 9h ~ aVminT. Pairs of Hawking 
particles are emitted close to the black hole horizon with 
a black body distribution at the Hawking temperature 
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The first equality is derived in Ref. [2|, the second one 
results from the explicit forms for v (9) and c (9) . 

Now, we can follow two different approaches to observe 
effects of Hawking radiation. In the first one we calcu- 
late the history of a final negative frequency pulse at 
small positive wave numbers leaving the black hole fol- 
lowing the standard techniques of l^, Ml- It is found 
that the final low wave number pulse scatters off the 
horizon and originates from two pulses, which conserve 
the frequency. These pulses consisting of high absolute 
wave numbers have such a low group velocity that they 
are moving rightwards in the lab frame, while being left- 
moving (upstream) in the comoving frame. The positive 
or negative wave number pulse has negative or positive 
frequencies. A different situation described in [IJ] ap- 
pears, when the two low wavelength solutions of the fre- 
quency condition ivq = vk ± D (k) , with the typical final 
frequency loq, lie outside the Brillouin zone. Then the 
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incoming pulses are instead rightmoving (downstream) 
and the positive or negative frequency pulse is located 
at high positive or negative wave numbers. This effect 
is analogous to Bloch oscillations. The particle produc- 
tion associated with Hawking radiation is found through 
the comparison of the so-called Klein-Gordon norm TV 
of the positive frequency initial pulse and the normalized 
final pulse. This initial pulse is found to reflect a thermal 
distribution at the Hawking temperature 13, | • 

We now return to the discrete ion chain. Our numerical 
results for the propagation of a final pulse backwards in 
time in the discrete system of phonons on an ion ring are 
presented in the following. We first introduce the quan- 
tities necessary for this analysis. If the phononic excita- 
tions are localized in the flat subsonic region, the excita- 
tions 66i(t) and their canonic conjugate S9i{t) can be ex- 
pressed as modes 59k(t) and 59k(t) with wave number k. 
Because of the finite system size only discrete wave num- 
bers appear and thus the dispersion relation is also dis- 
crete. The positive and negative frequency part of these 
excitations are defined by 269^ (t) — S9k{t) ± iS9k(t)/uJk 
and 269^ (t) — S9k{t) =F iuJkS9k(t). The analysis of the 
particle production requires us to use the Klein-Gordon 

69^*69^. We 
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norm for these modes J\f^ 



start from the final pulses 
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to test different frequency ranges. We calculate its dy- 
namics with the linear equations of motions following 
from the Hamiltonian (|4|). Our proposal does not suf- 
fer from the problem reported in J^] despite the steep 
dispersion at low wave numbers for long range Coulomb 
interactions. In our case, frequency conservation in the 
lab frame allows three pulses on the leftmoving branch 
of the dispersion relation exactly as in isl 13]. 

[13, [Til 
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Our results agree mainly with those of Refs 
shown in Fig. [51 We find the two previously mentioned 
high wave number pulses. In addition, we observe a weak 
rightmoving pulse with small negative wave numbers (see 
Fig. HI, whose presence agrees with approximate fre- 
quency conservation in the lab frame. This pulse could 
be a consequence of the finite initial and final excitation 
probability of ions outside the flat subsonic region. For 
small ion velocities the effect reported in [I^ analogous 
to Bloch oscillations is observed. Then the main incom- 
ing pulses are rightmoving. 

We further compare the Klein-Gordon norm of the pos- 
itive frequency initial pulse with the prediction for ther- 
mal radiation 13 li | . If only nearest-neighbor Coulomb 
interactions are considered, the relative difference be- 
tween these norms is lower than e — 0.01 for N — 1000 
ions. For the long range Coulomb interactions, a conser- 
vative estimate yields e < 0.2, which is in good agree- 
ment with thermal radiation. This conclusion is also 
confirmed by the analogous spectrally dissolved compar- 
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FIG. 2: Values of SOt{t) for propagation backwards in time 
starting from the final wave function in Eq. ((7| with s — 5. 
We use aVrmnT = 27r X 0.375, and iV = 1000. The fi- 
nal negative frequency pulses are depicted in green (solid 
line), the initial negative frequency pulses in blue (dash- 
dotted line), and the initial positive frequency pulses in red 
(dashed line). We depict the discrete norm distributions and 
dispersion relation as a continuous curve here, (a) Final 
{t = 0) and initial it = — 0.67r) norm distribution of 59k{t); 
eV47reo = 1.2591/(2iV) • mL^T-^. (b) Final {t = 0) and ini- 
tial (t — — 0.5ir) norm distribution of S9k{t) showing Bloch 
oscillation; e^/47reo = 2.0004/(2Af) • mL^T-'^. (c) Phononic 
group velocity in the flat subsonic region as a function of k 
for full interactions (blue dashed line) and nearest- neighbor 
interactions only (green solid line); as in (a). 



ison (see 'l3l|). Here the definition of a global Hawking 
temperature Th is difficult because of the nonlinear dis- 
persion relation at low wave numbers. We calculate Th 
from the discrete dispersion relation at low wave num- 
bers. The value of e can be explained with the inequality 
between the group and the phase velocity. Therefore, 
we have shown that Hawking radiation qualitatively per- 
sists even in a fundamentally discrete system with long 
range interactions and a logarithmically diverging group 
velocity at low wave numbers [see Fig. [5J:] . 

For an experimental proof for Hawking radiation on 
the ion ring, we propose to measure correlations between 
the different regions on the ion ring after the creation of 
a black hole following [5] . Starting from a thermal state 
of the excitations around homogeneously spaced ions at 
rest, the ions are brought into rotation and in a sub- 
sequent small time interval t a supersonic region is cre- 
ated. This is done by reducing the subsonic fluid velocity 
Vmin in a Gaussian way, while leaving the average rota- 
tion velocity constant. We are analyzing the normalized 
angle-angle correlations 

a,^{p.p,}-T/{hm). (8) 

In Fig. [3] we show the results at a fixed time starting 
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FIG. 3: Correlations Cij{t) at time t with real space lab frame 
positions. We consider A'^ = 1000 ions, aVminT = 2-k x 0.25, 
and T = 0.05T. (a) Only nearest-neighbor interactions 
are considered with e^/47reo = 1.127/(2iV) ■ rnL^T'^ and 
t — 0.5T. (b) Full Coulomb interactions are considered with 
eV47reo = 0.2453/(2Ar) ■ mL^T-'^ and t = 0.6T. 



from the ground state. Correlations between the super- 
and subsonic regions are created close to the black hole 
horizon and are moving away from it as expected. These 
correlations correspond to the pair creation of Hawking 
particles. We interpret their pure existence as a signature 
for Hawking radiation. Figure shows the simulation 
for interactions between neighboring ions only [l8| . Fig. 
[Sja for long range Coulomb interactions. The correlations 
behave similarly in both cases. We can tune and predict 
with good accuracy the angle and the propagation veloc- 
ity of the correlation features based on the group veloc- 
ities, i.e., the dispersion relation, and the ion velocities. 
The dashed lines in Fig. [3] show these predictions for the 
direction of the main correlation signal. 

We now discuss the parameter regime necessary in an 
experiment. Because the ion velocity must lie in the same 
order of magnitude as the phonon velocity in the pro- 
posed experimental setup, e'^/Aireo ~ 1/ {2N)mL^ /T'^ is 
required. This is satisfied for N = 1000 singly charged 
^Be ions with an average spacing of L/N = 10 fim and 
a rotation frequency of the ions of uJion — 2Tr x 5.5kHz. 
The Hawking temperature in this system is kBTn/h ~ 
5/T « 27r X A.ikHz. We have checked that the corre- 
lations of Fig. [3] remain present with initial tempera- 
tures 2 orders of magnitude above the Hawking temper- 
ature 0. At too high initial temperatures the emerging 
correlations are classical. The quantum nature of the 
correlations at low enough temperatures can be studied 
through the entanglement in the form of the logarith- 
mic negativity [l9j which is created between the inside 
and the outside of the black hole. We detect increasing 
entanglement in time between two regions adjacent to 
the horizon for initial temperatures up to 10 times the 
Hawking temperature. These parameters and this tem- 
perature should be reached in current experiments. We 
have checked that it is possible to control the equilibrium 



motion of the ions and to measure the correlations in the 
oscillations around them with sufficient accuracy. Note 
that it has been demonstrated long ago how to trap ions 
in quadrupole ring traps [l^ and measure their arrange- 
ment [lij . Thus, the proposed experiment will allow for 
the measurement of signatures of Hawking radiation for 
acoustic black holes. 

To conclude we have shown that it is possible to realize 
an acoustic black hole in a system of ions rotating on a 
ring. In this system thermally distributed Hawking radi- 
ation is expected, and we proposed a realistic experiment 
to measure signatures of Hawking radiation. We expect 
ring traps to be extremely useful for future quantum sim- 
ulations, especially of translationally invariant systems. 
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